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We argue that fermion-boson mapping techniques represent a natural tool for studying many- 
body supersymmetry in fermionic systems with pairing. In particular, using the generalized Dyson 
mapping of a many-level fermion superalgebra with the SU(2) type of pairing we investigate two 
kinds of supersymmetry connecting excitations in the systems with even and odd particle numbers: 
dynamical supersymmetry, which ensures a unified classification of states for both even and odd 
populations, and invariant supersymmetry with actual degeneracies of states within the same su- 
permultiplet. Dynamical supersymmetries based on the dynamical algebra of the type XJ(K/2Q) 
fC) t (where K and 20. denote the number of fermion-pair and single-fermion states, respectively) natu- 

rally arise in the bosonized description of the system. Conditions for invariant supersymmetry are 
' determined in a restricted case of bilinear supercharge operators. 
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I. INTRODUCTION 



^ ' Supersymmetry is commonly known as a hypothetical algebraic scheme in quantum field theory that unifies internal 
and space-time symmetries and predicts elementary particles of matter occurring in boson- fermion doublets [1]. As 
shown by Witten [2], the same scheme can be applied in nonrelativistic quantum mechanics, yielding a generalized 
dynamical algebra (superalgebra) which leads to an analogous doublet structure (except for the unique ground state) 
in the spectra of relevant quantum systems. In the simplest case, the SUSY Hamiltonian H = a) a + b^b is just the 
| anticommutator of supercharges Q = ba^ and = ab^ that change bosonic excitations (b') into fermionic ones {a') 
, and vice versa, so that the excited states \n a , rib) (where n a and rib, respectively, are numbers of fermionic and bosonic 
quanta) exhibit the characteristic two-fold degeneracy H\Q,rib) — H\\,rib — 1), coupling boson-like and fermion-like 
^3 , excitations. 

1 1 ■ Applications of SUSY quantum mechanics soon followed. Analytically solvable and isospectral sets of potentials 
were constructed on the basis of SUSY (see Refs. [3, 4] and references therein) and some of these potentials were 
found relevant for experimental spectroscopic data of certain atoms and ions. Thus an (approximate) manifestation 
of phenomenological supersymmetry was established in atomic physics [5, 6]. Methods based on the SUSY formalism 
were also developed in random matrix theory and applied to systems that exhibit signatures of quantum chaos [7] . 

In nuclear physics, the concept of supersymmetry found a natural application (see e.g. Refs. [8, 9, 10, 11]) in 
the framework of the Interacting Boson- Fermion Model [12]. The embedding of the IBFM dynamical algebra of 
the type Ub(-K^) ® Uf(2£!) (formed by K 2 bosonic generators b\bj and by (2fl) 2 fermionic generators a\ai, where 
i,j — 1, . . . K and k,l — 1, ...2fi enumerate single-boson and -fermion states, respectively) into the dynamical 
superalgebra \](K/2U) (with the mixed generators b\ai and atbj added) makes it possible to simultaneously describe 
low-lying spectra of doublets [8, 10] or quartets [11] of neighboring even and odd nuclei using a single energy formula 
with SUSY-based quantum numbers. (Note that in this paper direct products are used also for algebras, although 
rigorously we should speak about a direct sum of generators associated with the corresponding product of groups.) 

Dynamical supersymmetries are associated with chains T> Z> A\ D A2 D ... of possible algebraic decompositions 
of the dynamical superalgebra T> and refer to situations when the Hamiltonian H can be written solely in terms of 
the Casimir invariants C(Ai) of (super)algebras in one particular chain, thus implying complete integrability with a 
multitude of conservation laws [H, C{Ai)\ = 0. 

It should be stressed that nuclear supersymmetry, introduced by Iachello [8] already in 1980, was historically the 
first application of the SUSY ideas in nonrelativistic physics. Their most detailed verification in nature up to date — 
in the recent experimental work by Metz, Jolie et al. [13] — opened up further questions in the SUSY many-body 



'Electronic address: pavel.cejnar@mfT.cuni.cz 
^Electronic address: hbg@sun.ac.za 



2 



physics. Unfortunately, discussions of dynamical supcrsymmetry on the phenomenological nuclear-structure level, and 
its relation to the notion of SUSY quantum mechanics, have not always clarified the distinction between them, nor 
focused on the microscopic basis of nuclear SUSY. It is therefore not totally surprising to find that somewhat negative 
opinions such as the following one by 't Hooft [14] have been voiced: "At first sight, the fact that supersymmetric 
patterns were discovered in nuclear physics has little to do with the question of supersymmetry among elementary 
particles, but it may indicate that, as the spectrum of particles is getting more and more complex, some supersymmetric 
patterns might easily arise, even if there is no 'fundamental' reason for their existence." 

The dynamical supersymmetry described theoretically in Ref. [11], and found in experimental data [13], is based 
on an immediate decomposition 

U(Jf/2fi) D V B (K)®\] F (2Cl) D ... (1) 

of the dynamical supcralgebra into a product of the corresponding bosonic and fermionic algebras, namely on the 
decomposition U„(6/12) <g> U w (6/4) D Ubj,(6) ® Ubtt(6) £g> Uf^(12) ® Uf w (4) D ... in the observed case, where v 
and 7r, respectively, stand for neutron and proton realizations of the above (super)algebras. This implies that the 
achieved SUSY description does not really go beyond an application of the same IBFM Hamiltonian (with fixed 
interaction strengths) to a quartet of neighbouring even-even, even-odd, odd-even, and odd-odd nuclei. Indeed, since 
the dynamical-superalgebra irreps with more than one fermion in a given nucleus are highly excited, and thus do not 
mix with the low-energy spectrum, the SUSY Casimir invariants can be skipped as they only yield the same additive 
constants to level energies in all four nuclei considered [11, 13]. The remaining Casimir invariants (of the embedded 
algebras) are precisely those that appear in the standard IBFM description. 

As recently demonstrated in Ref. [15], any collective superalgebra of fermion-pair and singlc-fcrmion operators, 
linked to a microscopic perspective of the IBFM, can be naturally embedded into the phenomenological-type boson- 
fermionic dynamical superalgebra U(if/2f2). Here K and 2f2 now denote the number of pair and single- fermion 
operators, respectively, contained in the fermionic collective superalgebra. Moreover, the same analysis leads to 
the conclusion that all microscopically relevant dynamical supersymmetries must separately conserve the numbers 
of bosons and fermions, so that they are always of the form (1). Dynamical supersymmetry based on the IBFM is 
therefore always restricted to the structure (1), and should then probably only be referred to as dynamical boson- 
fcrmion symmetry. 

An interesting, though so far mostly hypothetical possibility was, nevertheless, considered by Jolos and von Brcntano 
[16]. It is based on the requirement that supercharges Qi, generating the odd sector of a boson-fcrmion superalgebra 
T C \J(K/2Si), together with generators of the even sector, commute with the Hamiltonian, and thus T forms a 
superalgebra of invariant supersymmetry of the system. This should be satisfied regardless of whether or not the 
system possesses any dynamical supersymmetry of the form (1). Clearly, the invariant supersymmetry represents 
principally the same kind of symmetry as the one introduced by the SUSY quantum mechanics, although the form of 
the Hamiltonian can be more general than simply the anticommutator of supercharges discussed above. It should be 
stressed that the bosons considered in this work are introduced on a purely phenomenological level. 

The distinction between the above two SUSY schemes is similar to the difference between invariant and dynamical 
symmetries in standard quantum mechanics: in the former case, all generators of the given symmetry algebra commute 
with the Hamiltonian, while in the latter this is satisfied only for Casimir operators of a certain chain of subalgebras 
(the last algebra in the chain thus represents an invariant supersymmetry). Dynamical supersymmetry, as presented 
in nuclear data, implies that all states in the neighbouring even and odd systems arc labelled by the same set of SUSY 
quantum numbers, but these states are not necessarily degenerate. The invariant supcrsymmetry, on the other hand, 
results in actual degeneracy (at least on a relative scale) of a certain subset of states, which has not been observed 
yet. Needless to say, that in field theory, as well, the SUSY scenario naturally assumes breaking of the supersymmetry 
between elementary bosons and fermions, as large mass differences for particles within the same supermultiplet must 
be accommodated. 

In this work, we extend the approach of Ref. [15] to study microscopic conditions for supersymmetric schemes in 
fermionic many-body systems. Our present discussion aims to clarify that while dynamical boson-fermion symmetry 
seems to be firmly established, and in the final instance based on the nuclear interaction in situations where some 
collective pairs are favoured (together with the utility of boson-fermion mappings, see Sec. II), the direct analogue 
in many-body systems of SUSY quantum mechanics (invariant supersymmetry), still requires clarification and may 
indeed be difficult to realise. This is explicitly demonstrated in Sec. III. Again the precise nature of the interaction 
on the fermionic level will of course turn out to be crucial for the appearance of invariant SUSY. 
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II. BOSONIZED FORMS OF THE PAIRING SYSTEM 
A. Pairing Hamiltonian 

In this section we present a bosonized form of a many-level fermionic system with the SU(2) type of pairing, a 
standard interaction to capture the essence of superconductivity in solids and nuclei [17]. We will deal with a set 
of single-particle states, enumerated by k — 1, . . . i7, and the respective time-reversal conjugate states, denoted by 
k = 1, . . . 0, occupied by an arbitrary number Np £ [0, 2f2] of fermions. The fermion creation and annihilation 
operators corresponding to individual states are a\, at and a k , a k . The Hamiltonian 

Hp = ^E k (^a\a k + a\aj^ -^G k ia\a\aiai (2) 
k k,i 

= Y.EkNk-Y.GuBlBi 

k k,l 

contains as free parameters single-particle energies E k (states k and k are degenerate due to the time-reversal invari- 
ance) and strengths of the pairing interaction G k i = G\ k . We can anticipate that the pairing interaction acts only 
between fermions in a certain subset V = {k ll k 2l ■ ■ ■ (with all fcj's mutually different, 5 < ft) of states — most 
likely states in some interval around the Fermi energy — so that G k i ^ only for k, I G V . 
The operators 

4 = 44 ( 3 ) 

and B k = {B\) + = a k a k create and annihilate, respectively, a pair of fermions (bifermion) on the shell k (thus there 
are i7 such pairs), and together with the /c-shell fermion number operator N k — a k a k + 4 a fc form the SU(2) algebra 

associated with each shell; [B' k , B k ] = N k — 1. Thus the Hamiltonian (2) leads to the dynamical algebra fc SU(2) fe . 
For special choices of E k and G k i, the Hamiltonian can be rewritten via a smaller number (K < 0) of pairs, 

k 

and = (C^) + (/i = 1, . . . , K), with the coefficients \nk satisfying the normalization constraint J2 k XvkXtk = 
The condition for the Hamiltonian reads as Hp = Hp, where 

H° = const + U„ v [Cl C v ] + VrvClC 




{Ufiv = U*,Vfj, v = V* are arbitrary interaction constants). Since creation and annihilation operators of any set 
of fermion pairs, together with their commutators, closes under commutation, the dynamical algebra generating the 
Hamiltonian in Eq. (5) can be identified with the algebra associated with the pairs (4) instead of those in Eq. (3). We 
call the C-bifermions in Eq. (4) collective. Note that they couple fermions into pairs in entangled single-particle states, 
in contrast to the noncollective B-bifermions (3) whose wave functions are separable (except, of course, unavoidable 
entanglement due to the antisymmetrization) . Note that Eq. (4) represents the most general ansatz for the collective 
pairs that conserves the form of the Hamiltonian (2). 

The possibility to rewrite the pairing Hamiltonian (2) via a collective algebra with a small number of C-bifcrmions 
is not very exceptional. In general, it requires to satisfy a set of equations for unknown coefficients x^k and for 
interaction parameters U^ Vl V^, cf. Eqs. (2) and (5). Dimensionality considerations suggest that some solutions 
should typically exist if K/Vt > yjl + l/tt + (E/fl) 2 — 1, which for S <C fi indeed leads to K -c ^- Let us stress 
here the importance of the requirement that also single-particle terms in Hp, not only the interactions, are expressed 
via the collective algebra, see Eq. (5). For instance, the constancy of interaction strengths in Eq. (2), G k i = G, 
would not imply that the simple SU(2) collective algebra generated by the single bifermion = J2kev 44/ v/ ^' ' s 
the dynamical algebra of the problem, unless all the single-particle states subject to interactions are degenerate, i.e., 
E k = E for k E V: 



Hi = ^ E k N k + EN V -GEC^C , 

k£V 



(6) 
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where N v = J^kev N k- 

Another example of the use of a collective algebra concerns a system of fcrmions on K single- j shells, j = 
3\,3i, ■ ■ ■ ,3 k, with a Hamiltonian expressed just via interactions of J = pairs on individual shells, — ( a j M a j M )^ : 

^=Ew-E^- ( 7 ) 

(E,j, are the j-shell energies, iV M the corresponding occupation-number operators, and V^„ the interaction strengths). 
Again, the Hamiltonian (7) is of the general form (2), with the conjugate states k and k corresponding to opposite 
projections +m M and — m M of the angular momentum for the same level jj,. 



B. Mapping via noncollective pairs 



The dynamical algebra of the fermion pairing problem can be translated into boson language with available boson 
mapping techniques [18], mostly explored for even fermion systems. If both even and odd particle numbers are to be 
described simultaneously, the algebra of fermion pairs must be extended to take into account also the odd degrees 
of freedom. This can most naturally be done by adding the single-fermion creation and annihilation operators to 
the given set of bifermions. The resulting collection of operators then forms a superalgebra with the algebra of pair 
operators as a subalgebra. 

As shown in Refs. [19, 20, 21], partial bosonization of this extended superalgebra can be achieved by a generalized 
Dyson mapping which, in the case of the superalgebra based on the above noncollective pairs (3), leads to 

B\ -> b{(l-n k -Af k ) , (8) 
B k >-+ b k , (9) 
N k i ► 2n k +M k (10) 



for the even sector, and 



, il-Af k -n k t t \-M k -n k 

bla k + al — — r a\a.\a k — rFT (H) 

k k k 1-A4 k k fc (2-A4)(l --A4) 



a k 



ak + a ' kh T^r k + al ^ akhk (2-N k )d-N k ) (12) 



for the odd sector [15, 21]. In the right-hand side images, bt (b k ) creates (annihilates) a boson of the fcth type and n k = 
b\.b k , while a\ (or a k ) create (or annihilate) ideal fermions in the states k (similarly for k), and M k = ot\a k + a k a k . 
These bosons and ideal fermions are kinematically independent, i.e. all boson operators commute with all fermion 
operators. The images of a| and a k can be deduced from Eqs. (11) and (12), respectively, after the k <-> k exchange 
of indices in the fermionic operators and inverting signs of the first term in Eq. (11) and of the second and third term 
in Eq. (12). 

The term 'ideal fermion' is used to distinguish the fermion-type particles resulting from the mapping (where they 
appear as necessary ingredients of any superalgebraic extension of boson mapping techniques) from the 'real' fcrmions 
in the original formulation of the problem. The Bogoliubov-Valatin like structure of Eqs. (11) and (12) suggests a 
physical interpretation of ideal fermions as generalized quasiparticles [22]. Anticipating discussions that follow below, 
we already point out here that possible divergences associated with the denominators in these formulas simply do 
not appear in the physical subspace. Furthermore the loss of symmetry between creation and annihilation operators 
under Hermitian conjugation (a typical feature of the Dyson type of mappings) is accounted for in the calculation of 
observables within the general Dyson framework, see Refs. [15, 18]. 

The mapped Hamiltonian reads as follows: 

H B = Y,E k (2n k +Af k )-^Gkib{bi(l + S kl -n k -Af k ) . (13) 
k k,i 

It is clear that this Hamiltonian conserves both the total number of bosons, Nb = J2k n k, an d °f ideal fermions, 
Afp = ^2 k J\f k , so that by considering a fixed total number of real fermions, N-p = ^2 k N k , we also fix the sum 
2N B +JVf to the given value iV F ; see Eq. (10). 

The Dyson mapping of states can be trivially deduced from Eqs. (8)-(12). Four real-fermion basis states in each 
A:-subspace yield the ideal boson-fermion images as follows: |0) ^ |0), a|.|0) ^ c4|0), c4|0) i— ► at|0), B' k \0) ^ b\.\0); 
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we denote real- and ideal-space vectors by angular and circled bras/kets, respectively, as in case of real and ideal 
vacua |0) and |0). All the states in the ideal Hilbert space that are not linear combinations of the images just given, 
are spurious. In particular, the components containing in any fc-shell more than one ideal particle (whether bosons 
or fermions, or both) are nonphysical. We see, therefore, that the present model of pairing allows us to write the 
projector to the physical subspace in the following explicit form: 

P p h = l[(Pn k =oP^ k =0+Pn k =oP^ k = l+Pn k = lP^ k =o) , (14) 
fe 

where P's on the right-hand side represent projectors onto the ideal subspaces with the given number of ideal particles 
(bosons or fermions) of the kth type. 

Because of the nonunitary character of the Dyson mapping, the boson-fermionic Hamiltonian in Eq. (13) (as well 
as images of other physical observables) is non-Hcrmitian with respect to the standard boson-fermion Fock space 
inner product. However, it is well known that Dyson mappings lead to so-called quasi-Hermitian structures which 
are consistent with standard quantum mechanics, and in particular guarantees real eigenvalues for observables (see 
Rcfs.[15, 23, 24]). Although explicit hermitization has been achieved in some particular cases [24], it seems that this 
procedure will generally introduce higher order many-body interaction terms into the Hamiltonian. Fortunately, in 
the present case, it can be easily checked that the Hamiltonian in Eq. (13) is already Hermitian within the physical 
subspace, i.e., it satisfies -ffe-Pph = H^Pph, where 

Hb = ^£ fe (2n fe +A4)-^G fe 4&* 

fe k,l 
= J2 E k^k+J2( 2E k- G kk)n k -J2 G kib{bi (15) 

fe fe fc^Z 

is manifestly Hermitian. We see that the Hamiltonian (15) contains no interaction terms — it is just a combination of 
bosonic and fermionic mean fields — although both boson-boson and boson-fermion interactions were present in the 
original mapped Hamiltonian (13). This implies that the interactions in Eq. (13) do not affect physical states. 

The Hamiltonian (15) has the same form as the one in Eq. (3) of Ref. [25]. Note that while in their treatment of 
pairing the authors of the cited work introduce boson-like particles, so-called cooperons, by modifying fundamental 
anticommutation relations of real fermions, our bosons result from the mapping of a conventional multifermionic 
superalgebra. In contrast to cooperons, the bosons in our case obey ordinary commutation relations, but for the 
limitations concerning the physical subspace they are in fact hard-core bosons, i.e., bosons with occupation numbers 
restricted to and 1. 

In spite of the single-particle form of the Hamiltonian (15), its diagonalization in terms of some new bosons 
d\ = 0i k b k and identification of the ground state for even systems with a condensate 

|cond)cx (5>/4J |0) (16) 

(where Nb — N-p/2 and rj k — (3i akl with Iq denoting the d-boson with minimum energy), arc not physically allowed 
operations. Indeed, such a procedure does not keep under control physicality of the transformed states; the condensate 
state has a finite spurious admixture as it contains terms with more than one boson on a given fe-shell (cf. Ref. [26]). 
The overlap of the state (16) with the physical subspace of the ideal space is given by the following expression: 

|(cond|F ph |cond)| 2 _ I Efc^fc^...^ Vk 1 Vk 2 ■ ■ ■ Vk NB ? 

) (.1 ') 



(cond|cond) ~ (£ fc \ m \ 2 f B A B ! 

which is less than unity for 1 < Nb < Np/2 and any nontrivial set of ?7's. At the same time, P p h|cond) is generally 
not an eigenstate of H B , nor of Hb- 

On the other hand, properties of the ground state of an even system can be estimated from the bosonic equivalent 
of the BCS approximation [17], using the image of the BCS state |BCS) = rife( u fe + v kBl)\0), namely 

|BCS) = l[(uk + v k bl)\0) = F ph exp^ (\nu k + —b{) |0) (18) 
fe k V Uk J 

with the normalization condition u 2 k + v\ = 1 for each k. The trial wave function (18) yields the following energy 
functional: 

(BCS|^ - 2AA B |BCS) - ^(2^ - G kk - 2\)v 2 k - £ G kl u k v k u m , (19) 

fe k^i 
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(with A denoting a Lagrange multiplier), the same expression as in the standard BCS approximation. 

We see that the bosonic BCS state (18) is just an ordinary Glauber coherent state projected onto the physical 
subspace (while the original BCS state is a generalized coherent state of the fermion dynamical group [27]). If the 
exponential on the right-hand side in Eq. (18) is projected onto a fixed number of particles (instead of the projection 
onto the physical space) one would get exactly the condensate state (16) with rjk = Vk/uk, but — as discussed above — 
this would not be a physically justified procedure. Instead, the use of P p h leads to the cutoff of spurious components, 
which also modifies the normalization factor in front of exp^ fc r)kb' k \0) in Eq. (18) with respect to the standard bosonic 
coherent state. 

C. Mapping via collective pairs 

If the pairing Hamiltonian can be expressed in the form (5), one can perform the fermion-boson mapping via the 
corresponding collective pairs. This is particularly simple in both the SU(2) special cases discussed above, see Eqs. (6) 
and (7): The C-bifermion from Eq. (6) yields a c-boson, i— > c^[l — (n + 7V)/S], while the S'-bifermions from Eq. (7) 
map onto s-bosons as follows: st [1 _ +J\f ll )/(j li + Note that the bifermion annihilation operators 

are in both cases mapped trivially just onto the corresponding boson annihilation operators, and the fermion number 
operators from the collective algebras onto analogous combinations as in Eq. (10). For instance, i— ► 2n M + A/^, 
where and Af^ stand for the boson and ideal- fermion number operators, respectively, of the /uth level in Eq. (7). 
It is clear that in this case the physical space is not restricted to the s-boson occupation numbers n M = and 1, but 
contains states with up to + | . The single- fermion operators, that in both of the discussed cases supplement the 
bifermion algebra to yield the corresponding collective superalgebra, are also mapped analogously as in the previous 
section, cf. expressions in Refs. [15, 21] 

The collective algebra of the bifermions (4) is not generally of the SU(2) type because different collective pairs, 
and CI, fi^y, may contain components with the same k. Nevertheless, using the formalism described in Refs. [20] 
and [15], the Dyson bosonic images can be easily constructed: 

X^kXvk X^kXvkXujkXirk c ljCir (20) 

v * k k ' 

h-* c„ , (21) 
[Cl,C v ] i ^ + ^] X^X^ ^ + 2 E E XMrt^X^Xrf 4 c t • (22) 

k lu.t; k 

Now ct and create and annihilate, respectively, the boson corresponding to the /xth collective bifermion. Note, 
however, that single-particle creation and annihilation operators generating the odd sector of the general collective 
superalgebra are more difficult to determine than in the SU(2)-based special cases [15]. 

It should be stressed that although the collective C-bifermions are just linear combinations of the noncoUective 
ones, = X^fcXjifc-Sfci the respective collective and noncoUective bosons are not connected by the same linear 
relation. Note that the two kinds of boson operators can be viewed to act formally in two different Hilbcrt spaces, 
so that their comparison requires to introduce an operator T which transforms the physical space of the collective 
mapping onto the physical space of noncoUective mapping. Eq. (21) and linearity of all mapping procedures trivially 
yield Tc^T^ 1 = = X^fcX^fc^fc! ^ means that the linear relation between C- and i?-bifermions is preserved for the 
respective boson annihilation operators. For the creation operators, however, the situation is much more complicated, 
as the application of T on the appropriate linear superposition of the mapping (20) results in a difficult selfconsistent 
relation containing combinations of operators = Tcj^T^ 1 , b\, bi, and also Afk = TAfkT^ 1 . This relation cannot be 

generally satisfied by any ansatz of the type = J2k X^fc^L where x'^k wou ld represent some unknown coefficients, 
and thus the correspondence between both kinds of mapping is not linear. 

III. SUPERS YMMETRY OF THE PAIRING HAMILTONIAN 

A. Possible SUSY schemes 

As explicitly demonstrated by the Hamiltonian in Eq. (15), the total number of particles is a conserved quantity in 
the system of bosons and ideal fermions obtained by the mapping (8)-(12). The same holds true also for any collective 
mapping, see Eqs. (20)-(22) (neither projection to the physical space, nor the Hermitization can spoil this property 
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[15]). Therefore, the dynamics of the mapped system does not rule out the use of the XJ(K/2Q) superalgebra as 
the spectrum generating (dynamical) superalgebra. In fact, the superalgebra of real fcrmions and bifermions, cither 
collective or noncollective, is mapped into (but generally not onto!) this superalgebra. 

For the noncollective mapping, for instance, K = Q, and the even sector of U(f2/2f2) is formed by generators 
b\bi, or k ai, a\aj, and their Hermitian conjugates (all together 5fl 2 operators), while the odd sector is generated by 
b\ai, b' k aj, and the conjugates (all together 4£1 2 ) operators. We see, however, that the odd sector is not used in the 
construction of the Hamiltonian, so that Ub(^) <£> Up(2£l) can equally well be chosen as the dynamical algebra. This 
means that the decomposition in Eq. (1) is applicable, which preselects possible dynamical symmetries of the system. 
Realizing that all these symmetries are given by standard decompositions of the Ub (K) £g> Uf (20) algebra into chains 
of embedded subalgebras, we skip here their explicit discussion. In general, conditions for any specific dynamical 
symmetry, given by a chain T> D Ai D Ai D ... of decomposition of the boson-fcrmion dynamical algebra T> 1 are 
obtained from the set of constraints [Hb, C(Ai)] = required to hold in the physical space. This naturally yields 
specific (for each symmetry) sets of equations for parameters Ek and Gki of the Hamiltonians (15) and (2). 

Although — as we just saw — the odd sector of the dynamical superalgebra U(K/2Q.) can be ruled out, it is also clear 
that the SUSY quantum number H = N-q + Np still classifies all eigenstates of the mapped many-body Hamiltonian. 
Conservation of the total particle number is a common feature of the generalized Dyson mapping. Let us stress that 
this simple conclusion was not at all clear in the early days of nuclear supersymmetry, especially in view of the fact that 
the real fermion number maps onto the number of ideal fcrmions plus twice the number of bosons, Np i— » Afp + 2Ab, 
which excludes the particle-number conservation law for real fermions from the explanation of the conservation of 
N. Indeed, as shown in Ref. [15], the conservation law for H (as well as the separate conservation of Nb and A/f) 
follows from a consistent choice of the real-fermion dynamical algebra, which must represent the interaction as well 
as single-particle terms of the fermion Hamiltonian. Therefore, it is only the possibility of fully algebraic formulation 
of the fermion pairing problem that automatically creates the U(A/2il)-based boson- fermion description, with all 
potential dynamical-supersymmetry chains included. 

Potential invariant supersymmetries of the pairing Hamiltonian (2) also depend on the the energies Ek and in- 
teraction strengths Gki- On the fermion level, all such supersymmetries are described in the framework of the most 
general superalgebra of multifermion operators, 

A v^X'\ =a l--- a kn a i ■ ■ ■ a L a K ■ ■ ■ a K> a K ■ ■ ■ a K> ' (23) 

where a given element belongs to the even or odd sector, respectively, according to whether the difference AAf = 
(n + h) — (n' + n') is even or odd. Operators of the even sector transform states with even and odd fermion numbers 
separately, while operators of the odd sector interconnect even and odd populations. The invariant SUSY takes 
place if the Hamiltonian (2) commutes with any of the operators (23) belonging to the odd sector, which implies the 
existence of a degenerate supermultiplct of states differing by a given odd number of AN-p fermions. 

In the Hilbert space of bosons and ideal fermions created by the mapping (8)-(12), similarly, one can consider the 
superalgebra of operators 

B fci,...fc„,fci...Wi, 7 l„ = ^fci,...fc„,fci,...fc„ . t _ b t b b (24) 

fej,...fe n / ,fc' 1 ,...fc^, fc' 1 ,...fc„/,fcj,...fe^, '1 'TO 'l ' m ' 

where A*}'"' is defined in the same way as A k k )' "' in Eq. (23), but with the real-fermion operators substituted by the 
corresponding ideal- fermion operators. As in the previous case, even and odd sectors of the superalgebra (24) are 
distinguished according to the difference AA/f between the numbers of creation and annihilation operators in A k )' ". 

Although the mapping in Sec. II B was performed only for a small subset of the above general fermionic superalgebra, 
the images of operators (23) in the physical space can be constructed as the corresponding products of the images 
in Eqs. (11) and (12) [19]. The superalgebra (23) is thus mapped to series of operators (24), i.e., onto the given 
superalgebra acting in the space of bosons and ideal fermions. Trivially, even and odd sectors of the superalgebra 
(23) are not mixed by the mapping — the image of an element with a given value of ANp has AA/f + 2AAb = ANp, 
where AAb = m — m! , see Eq. (24). Taking ANp = +1, for example, we get an image with terms corresponding to 
(AA/f, AA b ) = (+1, 0), (-1, +1), (+3, -1), (-3, +2), . . . 

This expansion would hardly be tractable in its general form, but fortunately, it can be considerably simplified if 
states with more than one ideal fermion are physically irrelevant. Indeed, such states correspond to real-fermionic 
states with broken pairs which are supposed to have relatively high excitation energies. Under this restriction, the 
invariant-SUSY supercharges connect states in doublets with A^p and Np + 1 fermions, and can be constructed just 
from the (AA/f, AA^e) = (+1, 0) and (—1, +1) terms. The first term applies in case that Np corresponding to the left 
system in the doublet is even, the second if Np is odd. To make contact with the standard notion of a supercharge, 
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we here only focus on the second term that changes a fermion into a boson (or vice versa). It reads as 

Qt = Y b k fki(b\ b h a L a n) {Qkiai + Qkiai) , (25) 
k,i 

where f k i are arbitrary functions of operators conserving separately the numbers of bosons and fermions, and qu and 
qui are some coefficients. 

In spite of the above restrictions, the commutator of the supercharge (25) with the boson-fermion Hamiltonian still 
remains completely unknown without a further specification of fki ■ This makes the general theoretical determination 
of invariant supersymmetries in fermionic systems a rather nontrivial task, where no analytic insight seems available 
so far. Nevertheless, the mapping procedure, as outlined in Sec. II, provides a natural framework for the classification 
and analysis of various special cases. In the next paragraph, for example, we determine conditions for the invariant 
SUSY generated by supercharges (25) with f ki = 1. 



B. Invariant SUSY with bilinear supercharges 

We will look for conditions that the mapped Hamiltonian from Sec. II B must satisfy to commute with the simplest 
bilinear supercharges of the form 

Q ] = bUlkiai + quaj) = ^ b\a k i , (26) 
k,i k.i 

where b\ and a; (or a{) are noncollective bosons and ideal fermions, respectively, from Eqs. (8) and (12). If the 
supercharge is assumed to commute with the time reversal, the coefficients q k i and q k i should be set equal (up to a 
phase factor), but this we do not generally require here. 

The assumption [H^, Q^] = (= [Q, Hq] + ) leads to the following equation to be valid for any k and I: 

(2E k — Ei — Gkk)qu — Y^ Gk m q m i — , (27) 

m^k 

and the same also with the change of g's to q's. For a fixed value of I, the condition for the existence of nontrivial 
solutions q k i and q k i reads as 



dct 



Gn — 2i?i + Ei G\2 G13 . . . \ 

G21 G22 — 2£/2 + Ei G23 

G31 G32 G33 — 2i?3 + Ei ... 

/ 



0. (28) 



This equation has rather clear physical interpretation: it says that diagonalization of the bosonic part of the single- 
particle Hamiltonian (15) gives one of the eigenvalues equal to the fermionic energy Ei. Mutual conversions of the 
corresponding boson and fermion thus do not change the total energy. One can apply Eq. (28) simultaneously to a 
certain subset of fermionic states, I G S, which is equivalent with taking q k i — qui = for I £ S. The supersymmetry 
then concerns only a part of the spectrum. 

An important special case of the above condition is obtained if one assumes supercharges with diagonal matrices of 
coefficients, q k i = qk$ki and q k l = Qk^kl- Eq. (27) then leads to the following simple constraints: G k i — for k ^ I e S 
and E k = G k k for k e S, the coefficients q k and q k being arbitrary for k e S and zero otherwise. We thus have 

h b = h svsy with 

H' SUSY ^J2 E ^+^ + J2 E ^ 2n k+^k)- ]T G kl b\h . (29) 

keS k<£S k,l<£S 

The interpretation of this solution is the same as above, with only the difference that the bosonic part of the Hamil- 
tonian (15) is already supposed to be diagonalized a priori, allowing for the diagonal constraint on the supercharge 
coefficients. By inserting the above constraints into the original Hamiltonian (2), one can verify the existence of 
degenerated supermultiplets with N k — 1 and 2 (for k e S) directly on the fermionic level. 

The Hamiltonian (29) is trivially invariant under the set of transformations b\a kl b\a k , a\a kl a k a k , ata k , b k b k 
(and their Hermitian conjugates) with k € S. These are generators of the X = <S> k ^s Ufc(l/2) superalgebra, which 
thus forms the invariant-SUSY superalgebra of the problem. This is so in spite of the fact that the dynamical 
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algebra of -ffg USY i n Eq. (29) is just an ordinary algebra T> C Ur(^) ® Uf(20) where seemingly no supersymmetry 
is involved. In view of the discussion above, the SUSY Hamiltonian obtained from the general solution of Eq. (28) 
in case of nondiagonal matrices of supercharge coefficients receives the same algebraic interpretation, but with the 
original bosons b\ substituted by the new ones, d\,, resulting from the diagonalization of the bosonic Hamiltonian. 
Note that the above superalgebraic scheme represents a direct generalization of the minimal SUSY scheme [2] , where 
the invariant-SUSY superalgebra sl(l/l) is formed only by operators Q\ Q, and H = {Q^,Q} [3, 4]. In the present 
case, in particular, the relevant part of the Hamiltonian is not just a supercharge anticommutator (or a sum of such 
terms), but a linear combination of Casimir invariants of X. 

Let us turn, at last, to the crucial question of physicality of the above supersymmetric transformations. It was so 
far completely ignored in this section but it is indeed very relevant since the supercharge operator (26) can transform 
some of the physical states out of the physical space given by the projector in Eq. (14). In particular, the supercharge 
operator acting within the physical space may produce a nonphysical state with (n k ,Af k ) = (1,1) or (2,0), this 
being not generally excluded — in case of nondiagonal matrices of supercharge coefficients by the action of a; or aj, 
associated with b\ in Eq. (26). Consequently, some of the supersymmetric transformations that, as found above, 
constitute an invariant supersymmetry of the Hamiltonian H' B may turn to be nonphysical. 

To examine this question, let us introduce a physical supercharge operator 

PphQtPph = Q f P P h (30) 

where 

& = Y h k & klPn k =oPN k =0 + Y h k & kkPn k =0 ■ (31) 

k^i k 

Similarly, we define 

Hb=Y Ek ^k + Y^ 2Ek - Gkk ) nk - Y G kibihPn k =vP^ k =o , (32) 

k k kjtl 

that represents another form of the Hamiltonian (13). The condition for the commutation of the supercharge with 
the Hamiltonian in the physical space then reads as [Hb, &]Pph — 0, which yields the following operator equality, 

Y( Ek ~ G kk)b\a kk + ^2(2E k -E x - Gkk)b\a k iPn k =oPN k =o 

k kjtl 

- Y G kib{auPn k =oP^ k =o - Y G kib\ai k P ril= QP^ l= Q 

k^l k^l 

+ Y G kib\ai m {P ni=1 - P ni= o)-P„ fe =o-PWi=o-PA4=o 

- Yl G kib f k b} r MmkPn m =oP^ m =o = , (33) 

k^l^m 

to be valid in the physical space. It is not difficult to see that Eq. (33) can only be satisfied by imposing the above 
SUSY conditions for supercharges with vanishing nondiagonal qu and q k i elements, see Eq. (29) and the text above. 
The Hamiltonian (29) is the only invariant-SUSY Hamiltonian with supermultiplet states contained entirely in the 
physical space that can be constructed through a bilinear supercharge in Eq. (26). The other solutions, see Eq. (28), 
are spurious, i.e., have no real counterpart on the level of the original fermion system. 

Let us finally note that the search for invariant supersymmetries would be much more involved for the mapping 
performed via collective bifermions from Eq. (4). In this case, we again have to check if the expression [Hb,Q }P p h 

vanishes, where = k cj Ji (q tJ , k a k + q^ k a k ) is a collective bilinear supercharge and P^ h a projector onto the 
physical subspacc obtained via the collective superalgebra. Unlike the noncollective case, boson-boson and boson- 
fermion interactions in the Hamiltonian Hb are generally relevant also in the physical subspace [15] and must be 
considered in the commutator. Moreover, the calculation is obscured by the fact that P^ h has a more complicated 
form in the general case than for the noncollective superalgebra. The above-discussed relation between collective and 
noncollective bosons, see the end of Sec. II C, indicates that the use of collective supercharges has a similar effect as 
the generalization of Eq. (26) to the form (25). 
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IV. CONCLUSIONS 

We have pointed out the need to distinguish between the well-established notion of dynamical supersymmetry in 
many-body physics (nuclear structure in particular), and invariant supersymmetry in many-body physics as a direct 
analogue of SUSY quantum mechanics. 

In contrast to phenomenological studies based on the IBFM, our search for the SUSY patterns is derived from 
purely fcrmionic level. Existing boson-fermion mapping techniques allow us to introduce the relevant boson and 
fermion degrees consistently, and facilitates our exploration of invariant supersymmetry in the simple case of pairing- 
like Hamiltonians. 

For these Hamiltonians, and with supercharges restricted to a bilinear form, it turns out that invariant supersym- 
metry is ruled out for pairing other than the simplest diagonal form, see Eq. (29). This is not in contradiction with 
the result of Ref.[16] where the boson interaction does not follow from a mapped fermion interaction. 

From a microscopic point of view, considering and accounting for the physical subspace, has been crucial for a 
complete analysis, as shown in detail in subsections II B and IIIB. 

Although we have not exhausted all possibilities that may result in invariant supersymmetry for generalized pairing 
in a fermion system, going beyond the present analysis seems to present a tough problem when pursued analytically. 

Finally, as argued in greater detail in Ref . [28] , and as is also clear from the present discussion, there is no reason to 
expect (as, e.g., in Ref. [29]) that transfer operators in the context of dynamical supersymmetry, should be restricted 
to the simplest supercharges appearing in invariant supersymmetry situations. In fact, for dynamical symmetry the 
appropriate transfer operators are mapped images of single fermion operators, as, e.g., in Eqs. (11) and (12). In 
our presentation these images may indeed be considered generalized supercharges, but in general they are not of the 
simple bilinear type. 

Acknowledgments 

P.C. thanks J. Dobes, J. Hosek, and D. Nosek for interesting discussions. This work was supported by the S.A. 
National Research Foundation under Grant Nos. GUN 2047181 and GUN 2044653 and by the Grant Agency of Czech 
Republic under Grant No. 202/02/0939. 



[1] S. Weinberg, The Quantum Theory of Fields, Vol. Ill, Supersymmetry (Cambridge University Press, Cambridge, UK, 
2000). 

[2] E. Witten, Nucl. Phys. B188, 513 (1981). 

[3] F. Cooper, A. Khare, and U. Sukhatme, Phys. Rep. 251, 267 (1995). 

[4] F. Cooper, A. Khare, and U. Sukhatme, Supersymmetry in Quantum Mechanics (World Scientific, Singapore, 2001). 
[5] V.A. Kostelecky and M.M. Nieto, Phys. Rev. Lett. 53, 2285 (1984); Phys. Rev. A 32, 1293 (1985). 
[6] V.A. Kostelecky, M.M. Nieto, and D.R. Traux, Phys. Rev. D 32, 2627 (1985). 

[7] H.-J. Stockmann, Quantum Chaos. An Introduction, (Cambridge University Press, Cambridge, UK, 1999). 
[8] F. Iachello, Phys. Rev. Lett. 44, 772 (1980). 

[9] A.B. Balantekin, I. Bars, and F. Iachello, Nucl. Phys. A370, 284 (1981). 
[10] A.B. Balantekin, I. Bars, R. Bijker, and F. Iachello, Phys. Rev. C 27, 1761 (1983). 
[11] P. Van Isacker, J. Jolie, K. Heyde, and A. Frank, Phys. Rev. Lett. 54, 653 (1985). 

[12] F. Iachello and P. Van Isacker, The Interacting Boson-Fermion Model (Cambridge University Press, Cambridge, UK, 
1991). 

[13] A. Metz, J. Jolie, G. Craw, R. Hertenberger, J. Groger, C. Giinther, N. Warr, and Y. Eisermann, Phys. Rev. Lett. 83, 
1542 (1999). 

[14] G. 't Hooft, The Holographic Principle, Opening Lecture at the 1999 Erice-Chalonge School, hep-th/003004. 
[15] P. Cejnar and H.B. Geyer, Phys. Rev. C 65, 044313 (2002). 

[16] R.V. Jolos and P. von Brentano, Phys. Rev. 60, 064318 (1999); 62, 034310 (2000); 63, 024304 (2001). 

[17] J. Bardeen, L.N. Cooper, and J.R. Schrieffer, Phys. Rev. 106, 162 (1957); 108, 1175 (1957). 

[18] A. Klein and E.R. Marshalek, Rev. Mod. Phys. 63, 375 (1991). 

[19] J. Dobaczewski, F.G. Scholtz, and H.B. Geyer, Phys. Rev. C 48, 2313 (1993). 

[20] P. Navratil, H.B. Geyer, and J. Dobaczewski, Phys. Rev. C 52, 1394 (1995). 

[21] P. Navratil, H.B. Geyer, and J. Dobaczewski, Nucl. Phys. A 607, 23 (1996). 

[22] T. Suzuki and K. Matsuyanagi, Prog. Theor. Phys. 56, 1156 (1976). 

[23] F.G. Scholtz, H.B. Geyer, and F.J.W. Hahne, Ann. Phys (NY) 213, 74 (1992). 

[24] G.K. Kim and CM. Vincent, Phys. Rev. C 35, 1517 (1987). 



11 



[25] E. Celeghini and M. Rasetti, Phys. Lett. 283A, 382 (2001). 
[26] J. Dobes, Phys. Lett. 222B, 315 (1989). 

[27] W.-M. Zhang, D.H. Feng, and R. Gilmore, Rev. Mod. Phys. 62, 867 (1990). 

[28] H.B. Geyer and P. Cejnar, in Proceedings of the Eleventh International Symposium on Capture Gamma-Ray Spectroscopy 

and Related Topics, ed. J. Kvasil, P. Cejnar, and M. Krticka (World Scientific, Singapore, 2003), p. 43. 
[29] J. Barca, R. Bijker, A. Frank, and G. Loyola, Phys. Rev. C 64, 064313 (2001). 



